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Abstract. A class of periodic soliton cellular automata is introduced associated with 
crystals of non-exceptional quantum affine algebras. Based on the Bethe ansatz at q = 0, 
we propose explicit formulae for the dynamical period and the size of certain orbits under 
the time evolution in the A„ case. 



1. Introduction 



The box-ball system |T5HT] is a soliton cellular automaton on a one dimensional lattice. 
It is an ultradiscrete integrable system TTMS that exhibits factorized scattering, and has 
been studied from a variety of aspects. Among them an efficient viewpoint is a solvable 
vertex model in statistical mechanics Q5] at q = 0, where the time evolution of the box-ball 
system is identified with the action of a transfer matrix. It has led to a direct formulation 
HHIKTTl IFOYj by the crystal base theory a theory of quantum group at q — K , 
and generalizations associated with quantum affine algebras HKT1, HKO TY] . For some 
latest developments along this line, see [IKOl IKO Y| . These studies are based on the idea 
of commuting transfer matrices Q5]. As a method of analyzing solvable lattice models, it is 
complementary to the most efficient technique known as the Bethe ansatz |Be| ; therefore it 
is natural to seek its application to the box-ball system and its generalizations. 

The aim of this paper is to extend the box-ball system to periodic versions and launch 
a Bethe ansatz approach to them. For non-exceptional affine Lie algebra 0„, we construct 
a periodic ultradiscrete dynamical system that tends to the q„ automaton HKTl in an 
infinite lattice limit. Here is an example of the time evolution pattern for g rl = A^p: 



t = 
t = 1 
t = 2 
t = 3 
t = 4 
t = 5 
t = 6 
t = 7 
t = 8 
t = 9 
t =10 
t =11 
t =12 



112 13 2 
3 2 12 11 

113 12 2 
2 2 13 11 

1 1 2 2 3 1 

2 1112 3 
13 2 112 

2 113 2 1 
12 2 113 

3 112 2 1 
2 3 1112 
1 2 3 2 1 1 
112 13 2 



Regarding the letter 1 as background, one observes two solitons proceeding cyclically to the 
right with velocity = amplitude equal to 2 and 1. They repeat collisions (or overtaking) under 
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which the reactions 32 x 2 — > 3 x 22 and 22 x 3 — > 2 x 32 take place. Behind such dynamics 
there underlines a solvable vertex model at q = 0, where only some selected configurations 
have non-zero Boltzmann weights and the transfer matrix yields a deterministic evolution of 
the spins on one row to another. For instance, the transition from t — to t — 2 states has 
been determined from the configuration in figure 1 on a two-dimensional square lattice: 

112 13 2 

123 — | — 112 — | — 111 — | — 112 — | — 111 — | — 113 — | — 123 

3 2 12 11 

111 — | — 113 — | — 123 — | — 112 — | — 122 — | — 112 — | — 111 

113 12 2 

Figure 1 . A vertex configuration 

This is a configuration of the fusion U q (A^) vertex model that survives at q = Q. In 
the terminology of quantum inverse scattering method STFj, the quantum space on vertical 
lines carries the fundamental representation (1, 2 or 3) of A2 = sl^ and the auxiliary space on 
horizontal lines does the three-fold symmetric tensor representation (111, 112, . . . , 333). The 
automaton states live on the vertical lines. The dynamics is governed by combinatorial R, 
which is the quantum R matrix at q — specified by local configurations round a vertex. The 
states on horizontal lines are so chosen that they become equal at the both ends reflecting 
the periodic boundary condition. 

In this paper we introduce analogous periodic automata for any non-exceptional affinc 
Lie algebra Q n based on the factorization of the combinatorial R H KT2j . They may be 
viewed as the system of particles that undergo pair creation and annihilation though the 
collisions. Moreover we exploit how the Bethe ansatz at q = |KN| yields the dynamical 
period and size of certain orbits. For instance in the above time evolution pattern, t = 
and t = 12 states are identical, hence the dynamical period is 12. We propose the general 
formula l|16(l for the dynamical period in An case, which indeed predicts 12 in the above 
example. It is expressed as a least common multiple of the rational numbers arising from 
Bethe eigenvalues at q = 0. 

In KN , the Bethe equation is linearized into the string centre equation and an explicit 
character formula (|23|l has been established by counting off-diagonal solutions to the string 
centre equation. It is a version of fermionic formulae and is called the combinatorial com- 
pleteness of the Bethe ansatz at q = 0. In (|24Jl we relate each summand (1221 in the character 
formula to the size of a certain orbit under the time evolution. Such a result will be useful 
to study the entropy of the automata. 

The formulae for the dynamical period (|16l) and the orbit size l|24[) are novel applications 
of the Bethe ansatz to ultradiscrete integrable systems. Upon identification of strings in 
the Bethe ansatz with solitons in the automata, they reproduce the expressions in |YYT| 
for A± with I = 00. In our approach, we also use the combinatorial Bethe ansatz at 
q = 1 KKR, KR , namely, rigged configurations and their bijective correspondence with 
automaton highest states. In this terminology, it is the configuration that plays the role of 
the conserved quantity, which is an analogous feature to the infinite system [KOTYj . It is 
an interesting problem to synthesize the combinatorial Bethe ansatze at q = 1 and q = 0, 
which will provide a unified perspective on the automata on the infinite and the periodic 
lattices. 
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The paper is arranged as follows. In section a periodic Q n automaton is introduced. 
It tends to that in |HKT1| in an infinite lattice limit and includes the one in |YTL IMIT] as 
a Q n = An" 1 case. In section[31 the Bethe eigenvalues are investigated at q — 0. In section^] 
the dynamical period of the periodic automata is related to the Bethe eigenvalue studied in 
section [3J In section sizes of certain orbits are related to the character formula in |KN| . 
In the last two sections conjectures are presented with compelling experimental data. The 
states treated there are time evolutions of the highest ones. The classes of time evolutions 
being considered in sections 0] and are different. In fact, the latter is wider containing 
the former; therefore, the 'period' in section 21 is a notion different from the 'size of orbit' 
in section [5] A more unified framework including a treatment of non-highest states will be 
presented elsewhere. The last table in section0]is a preliminary report on For standard 

notations and facts in the crystal base theory, we refer to [Kj IKKML IKMN| . 

2. Periodic g n automaton 

Let Uq(g n ) be the quantum affine algebra associated with non-exceptional g n = An\ 

A fl-n A 2n> B n \ D { n ] and D^l v D enote by B t the crystal of the Wold symmetric 

fusion of the vector representation of U q {Q n ) jKKMj . We are going to introduce a dynamical 
system on the finite tensor product B := B{ x ® Bi 2 g) ••• Bi L . An element of B will be 
called a path. The representative time evolution is given by 



(1) 



Too — GFiSi d ■ ■ ■ Si 2 Si 1 . 



Here Si is the Weyl group operator [K] and ob = cr (8 • • ■ <3> o with each a acting on the 
components B\ i individually according to Table For example for the element 11245 € B$ 



of represented by the semistandard tableau, one has ct(11245) = 13455. Sec HKT2 
section 2 for the notation in the other algebras. We call the dynamical system on B with 



the time evolution JJ) the periodic Q n automaton. In case B is of the form B 
be called the basic periodic g n automaton. 



Bf ^, it will 



Table 1. The data a, d and ik 



A (2) 
A 2n-l 
A (2) 

B m 
r w 

D (2) 



d 


id, 


W 


n 




2, 3, ... , n—1, n, 




2n-l 


0,2, 


. . , n—1, n, n — 1, . . . , 2, 




2n 


1,2,. 


. ,n — l,n,n — 1, ... ,2, 1,0 


W{BC n -i) 


2n-l 


0,2, 


. . , n— 1, n, n— 1, . . . , 2, 


W(BC n -i) 


2n 


1,2,. 


. ,n— l,n,n — 1, ... ,2, 1,0 


W{BC n -i) 


2n-2 


0,2,... 


2 («-!,«) n _2,...,2,0 

' l n,n— 1 J ' ' ' ' 


W(Dn-l) 


2n 


1,2,. 


. ,n— l,n,n— 1, ... ,2, 1,0 


W{BC n -i) 



a i — ► a — 1 

1 <-> T 

id 
1 <-> T 

id 

1 <-> 1, n <-» n 

id 



Let us illustrate Q along a g n — A^ 1 example. The time evolution = crs^S'o of 



the t = path 112132 G Bf 6 into 321211 at t = 1 in section[T]is computed as 
(2) 



11 2132 Sl 3 2 3 32 Sl 3 2 3 2 2 K 3 2 1 2 1 1. 



— signature 
2— signaurc 



(+ 



(+-) 



(+-) + + 



For the first three paths, we have exhibited the 0-signature and 2-signature. In general, 
the i-signature of an element a in the A% crystal B\ = {1,2,3} is the symbol + if a = i, 
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— if a = i + 1 mod 3 and none otherwise. From the array of «-signatures, one eliminates 

a 

the pair H — (not — h) successively to finally reach the pattern — • • • — !-■■■ + called the 
reduced i-signature. Then the action of Si is unambiguously defined as the interchange 

a (3 Pa 

— ■ ■ ■ — !-•••+ i— > — ••■ — !-•••+ on the reduced i-signature. In @, we have shown the 
elimination of the H — pairs by parentheses. By the very same rule, the Weyl group operators 
in general g„ and B = Bi 1 ® • • • ® Bi L can be computed using the necessary data on B\ in 
KKM. 

One may wonder the relation between the two derivations of the time evolution 112132(i = 
0) i — > 32121 l(t = 1), one as @ and the other as in Figure ^ m section ^ Let us clar- 
ify it by explaining the origin of (Q. Recall that the automata in the infinite system 
HKT1, HHIKTTj IFOYL HKOTY have the set of states • • • <g) B h <g) B h+1 <g> • • • with 
the boundary condition that the sufficiently distant local states are the highest element 
ui i = e Bi t . The commuting family of time evolutions T\ (I 6 Z>i) is induced by the 
relation 

Bi <g> (• • • <g> B h ® B h+1 ® • • • ) ~ (• • • <g> B u <g> B h+1 ® • • • ) ® B t 
^ ui®p ~ Ti{p)®ui 

under the isomorphism of crystals. It was proved in HKT2 that 2] with sufficiently large 
I is factorized as (Q, where all the Si actually act as e°°. In this sense is a natural 
analogue of the in the infinite system, which corresponds to the limit of the periodic g„ 
automaton when the system size L grows to infinity under the above-mentioned boundary 
condition. The product is a translation in the extended affine Weyl group. The indices 
ik in Tableware equal to ik+j in HKT2 for some j. They have been chosen so that the 
tableau letter representing the background or 'empty box' to becomes 1. 

Let us comment on the analogue of the time evolution 7] with finite / on our periodic g„ 
automaton. A natural idea is to define it by an analogue of the relation 10 as vi®p ~ p' (givi, 
where vi € Bi is not necessarily the highest element ui in general. If such a vi exists and p' is 
unique even when vi is not unique, we set Ti(p) — p' and say that Ti(p) exists. Ti(p) does not 
always exist. For instance in An^ case, v\ does not exist for L = I = l,p = 12 £ B = B2, and 
p' is not unique for L = 2, 1 = l,p = 12 (g) 12 € B2 (8 B 2 - See section |S] for more arguments. 
On the other hand for / sufficiently large, we expect that Ti{p) exists. In fact the following 
assertion is valid. 

Theorem. Let g„ = An (hence d = n). Pick any element p G B such that 

(4) Vik (S ik _ 1 ---S il (p))<e ik (S ik _ 1 ---S il (p)) forl<fc<d. 

Set vi — (xx, . . . , x„+i). Here the number Xi € Z>o (i G Z n+ i) of the letter i in the semistan- 
dard tableau on length / row is determined by Xi k = ipi k (Si k l ■ ■ ■ Si t (p)) for 1 < k < n and 
x\ + • ■ • + x n+ \ = /, which is possible for / large. Then for sufficiently large /, the relation 

(5) vi <g)p~ Too(p) ®vi 

holds under the isomorphism of crystals Bi ® B ~ B ® Bi with given by Q • 

The condition stated in an intrinsic manner, is actually a simple postulate that 
among {1, . . . , n+1}, the letter 1 should be no less than any other ones in the semistandard 
tableaux consisting of p. On the paths in FigureUJ (112132, 321211, 113122 e Bf 6 ), one has 
Tn = T^. A similar theorem is valid also for DnK 
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133 


■A- 


1 


3 


,4. 


12- 


4 


124' 


3 


4- 


1 


3 


14. 


3 


134' 


2- 


3 


•4- 


1 • 


34- 


1 


134' 


1 • 


■ 2 ■ 


■ 2 ■ 


,4. 


13- 


4 



The time evolution Too Q commutes with several operators acting on _B, which form the 
symmetry of (Too flow of) our periodic g„ automaton. By using Q and Sicrs = <tbS -i(,) 
(see |HKT2| ). it is easy to check that T^Si — S'iTx, for i ^ 0, 1. The symmetry operators 
or 'Backhand transformations' {Si | 2 < i < n} form a classical Weyl group listed in the 
rightmost column of Tabled This is a smaller symmetry compared with the C^j(0„_i)- 
invariance in the case of the infinite system HKOTY] . 

Here is an example of the time evolution (downward) in the periodic A3 automaton 
with B = B 3 (gi Bi <g Bi (g) Bi (g B\ <g> B 2 <g> -Bi . At each time step, the paths connected by 
the Weyl group actions S2 and S3 are shown, forming commutative diagrams. 

h4 123-4- 1-2-4- 12-4 & 123 • 4 • 1 ■ 2 • 3 ■ 12 • 3 
124 • 3 ■ 4 ■ 1 • 2 ■ 14 • 2 123 • 3 • 4 • 1 • 2 • 13 • 2 
124 • 2 ■ 3 ■ 4 • 1 ■ 24 • 1 123 ■ 2 • 3 • 4 • 1 • 23 • 1 
124 • 1 ■ 2 ■ 3 • 4 ■ 12 • 4 123 ■ 1 • 2 ■ 3 ■ 4 • 12 • 3 

A similar example from the basic periodic D± automaton with B = Bf 12 . 

222211114211 h4 232211114211 & 242211113211 

111122221421 111123221421 111124221321 

211111112241 211111112341 211111112431 

223421111131 222421111131 222321111141 

121112342113 121112242113 121112232114 

413211111232 413211111222 314211111222 

331134321111 321134321111 321133421111 

133111113432 132111113432 132111113342 

Let us remark on another family of maps on B, which may also be regarded as time 
evolutions. For it is a dual of JJJ (cf. |KNY| ). Consider the maps Ti, . . . ,T L ; B — ► B 
defined by 

% = Ri-l i ■ ■ ■ R23Rl2PiRh-l L ' ' ' Ri+1 HQ.Ri H-l , 

(6) Pi : B yi <g> B h -» B h <g B Vl 

p®b 1 ^ b® p. 

Here £+i is the combinatorial R that exchanges the fc-th and (fc+l)-th component from the 

left, and B Vl = Bi x % ■ ■ v • • (g) B; L is the B without the component Bi t . It is an observation 
going back to \7\ that the Yang-Baxter equation and the inversion relation of R lead to the 
commuting family %Tj = TjTi- Note that % — %+% when £j = ij+i. 



3. Bethe eigenvalues at q = 

In this section we exclusively consider the simply laced g„ . Eigenvalues of row transfer 
matrices in trigonometric vertex models are given by the analytic Bethe ansatz IKS| . In 
the present case, the relevant quantity is the top term of Aj (it) ((2.12) in KS modified 
with a parameter h to fit the notation here): 

{ ' Qi(u + lh) 

at the shift (or Hamiltonian) point u = 0. Here Q\{u) — J\ k sinh.7r(u — y/— lu^ ), where 
{u^} are to satisfy the Bethe equation eq.(2.1) in |KN| . For the string solution ( KN 
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Definition 2.3), 10 with u = tends to 

(8) A i :=n(-^ ) ) min °'' i) 



in the limit q — exp(— 27th) — > 0. Here is the center of the a-th string having color a 

(a 
1 

j € Z>i and 1 < a < m^. At q — the Bethe equation becomes the string centre equation 
([EM (2-36)): 

(&) 

o) n n^?)^^-^ 1 )^^' 4 " 1 ' 

{b,k)£H 13=1 

where H := {(a,j) | 1 < a < n, j € Z>i,TOj- > 0} (denoted by £T in |KNj ). A a j a> bkp and 
are defined by 

(10) A a j a ^kp = SabSjkSapipf 1 + m^) + C a b min(j, fc) - 6 a b8jk, 

(11) p^ a) = ^ min(j, fe)^ o) - C a6 min(j, k)mf\ 

k>l (b,k)eH 

where (C a i,)i<a,b<n is the Cartan matrix of the classical part of g n . The integer v)f is 
the number of the Kirillov-Reshetikhin modules wjf^ contained in the quantum space on 
which the transfer matrices act. In our case, the crystal of the quantum space is taken as 
B = Bi x ® ■■■ ® Bi L in section [21 hence i/^ = S a i(Ski 1 + ■ • • + §kl L )- To avoid a notational 
complexity we temporally abbreviate the triple indices aja to j, bk(3 to k and accordingly 
z kp to Zk etc - Then (JSJ) reads 

(12) A, = n(-**r. 

fe 

where pu is actually dependent on I, and given by pk = Sbi min(&;,/) for fc corresponding to 
bk/3. The string centre equation © is written as 

(is) n(-zk) A ^ = (-1)^- 

for some integer Sj. Note that Aj t k — Ak.j. Suppose that the q = eigenvalue 1121) satisfies 
A^' = ±1 for generic solutions to the string centre equation i|13|) . It means that there exist 
integers rj such that J^j r jAj,k = Vipk, or equivalently 

, x det A\j] 

where A[j] denotes the matrix A = (Aj k) with its j-th column replaced by (pi, P2, ■ ■ •)■ I n 
view of the condition Wj 6 Z, the minimum integer value allowed for Vi is 

(15) ~ ' " <H 1 



and length j. Denote by rrij the number of such strings. The product in |JSJ is taken over 



where LCM stands for the least common multiple and L)' k means the union over those k such 
that A[k] 0. The determinants here can be simplified by elementary transformations (cf. 
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|KN| (3.9)). The result is expressed in terms of determinants of matrices with indices in H: 

where the matrix F = (F a j tb k)( a ,j),{b,k)£H is defined by 

(17) Faj,bk = SabSjkPj + C ab min(j, k)m%' . 

The matrix F[b, k] is obtained from F by replacing its (b, k)-th column as 

as) F [b ,k Ucm j^ { r\*%n 

^ai mm(j,() (c,m) = [b,k). 

The union in (|16|l is taken over those (b, k) such that detF[6, fc] ^ 0. See also the remark 
before Conjecture 1 in section 

The LCM in (|16() can further be simplified when g„ = A± and vf' = L5j\. We write 



in 



(i) Ji) 



. ,Pj ,F[l,k] just as nij,Pj,F[k] and parameterize the set H = {j 6 Z>i | mj > 0} as 
if = {(0 <) Ji < • • ■ < J s }. The matrix F[/c] is obtained by replacing the k-th column of F 
by *(min( Ji, I), min( J2, Z), . . .). A direct calculation leads to 

(19) detF=pj pj 1 ■■■pj,_ 1 , 

(20) det F[k + 1] - det F[fc] = PM ' ' ' PJ ^ " ih) , 0<k<s-l, 

Pik+lPik 

where we have set if. — min(Jk,l), Jo = 0, Iq = 0, po = L and F[0] = 0. (ik here is not 
related to those in Tabled) Substituting (fTTTfl into lfTH|) and using the elementary property 
of LCM, we find 

7>, = LCM(l, d6t ^ d6tF d6t ^ 



det F[l] ' det F[2] ' " ' ' det F[s] . 

_ / detF detF detF 

( 21 ) ~ V ' det F[l] ' det F[2] - det F[l] ' " ' ' det F[t + 1] - det 



= LCM 1, y ' 



where < t < s — 1 is the maximum integer such that i t +i > it- 

4. Dynamical period 

In this section we shall exclusively consider case although the parallel results are 
expected for D„ . When I — > 00, one puts ik = Jk and t = s — 1 in the formula (|21[1 . 
Eventually the resulting expression coincides with eq.(4.24) in |YYT| . which gives the period 
of generic paths in the periodic box-ball system containing mj solitons of length j. Here by 
generic is meant the absence of an 'effective translational symmetry' |YYT| . In the present 
framework it corresponds to the time evolution T; =oc of the basic periodic A^ 1 automaton, 
i.e., Q n = A^\v^ = LSji. 

To generalize such a connection, we invoke the combinatorial version of the Bethe ansatz 
explored in KKR , KR . Given a highest path, namely an element p G B such that iip = 
for l<i<n, one can bijectively attach the data (rrv- a \ r^ a ')2—i called rigged configuration. 
Here to'"' = (m^) is a Young diagram involving rows of length j, and — (r^) 
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stands for an array of partitions attached to each 'cliff' of 

m («). | m (o)| i s 

equal to the 

number of letters a+1, a+2, . . . , n+1 contained in the corresponding path p. The separate 



data (m 



(i) 



, m^™)) and 



"(")) are called configuration and rigging, respectively. 



They obey a special selection rule originating in the string hypothesis. Namely, the p 
defined by must be nonnegative and the maximum part of the partition is not 



greater than j^ a ' for any (a,j) e H. It is known that detF > ([KNj Lemma 3.7) for any 
configuration. Time evolutions of a highest path is not highest in general. Let Ph{m) C B 



be the set of highest paths whose configuration is m = (m 



Conjecture 1. For a highest path p € Ph(m), suppose that T^{p) exists for any k £ Z>i. 
Then the dynamical period of p (minimum positive integer k such that T, (p) — p) is equal 
to Vi (|16|l generically, and its divisor otherwise. 

Naturally we expect A^ 1 = 1, which can indeed be verified for . Conjecture 1 implies 
that the generic period is a function of the configuration only and does not depend on the 
rigging. We have abruptly combined the Bethe ansatz results in two different regimes. The 
first one in section [3] is at q = KN , whereas the second one explained here is relevant 
to q = 1 KKffJ IKRJ . Conjecture 1 has been confirmed for all the highest paths in Bf L 
with L < 9 and all the sl n <4 highest paths in Bi ± ® • • • ® Bi L with li + ■ ■ ■ + II < 7. It was 
observed that non-generic cases are pretty few and Tf(p) (k> 1) exists for any highest p in 
c&seB = Bf L . 

Let us present a few examples of Conjecture 1. To save the space 12 ® 224 is written as 
12 • 224 etc., and furthermore, • is totally dropped for the basic periodic automata. In each 
table, the period under 7] with maximum I is equal to that under Too. The last table is a 
preliminary report on the case. 



4 1} 



path 
I 



1211122122111221122, configuration = ((32211)) 



LCM of 



19, 19, 

57 i^i 

' 22 ' 

171 — 

5 22 ' 



19 

171 
22 
513 
193 



= period 



19 
171 
513 



A[ 1} path = 11 ■ 1112 • 2 • 112 • 122 ■ 2 • 2 • 1, configuration = ((43)) 



I 


LCM of 


= period 


1 


1 21 


18 


54 


2 


1 21 


9 


27 


3 


1 2 

x ' 2 ' 


6 


18 


4 


1, 9, 


3 


9 



A ( 2 ] path = 121121213322111133211, 



LCM of 



configuration = ((431 11), (4)) 
= period 



1, 21, 



21, 



21. 



21 



2 822 822 411 411 



21 

822 
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a£ ] path = 1 • 12 • 3 • 114 • 1 • 2 • 22, 



Z 



1. 



LCM of 



1 29 

•"■> 5 ' 



29 29 



29 
12 ' 



58 
9 ' 



29 
4 



i 58 58 116 116 
J-i 7 j 13 j !7 ; 17 



5s 
9 



configuration = ((3111), (11), (1)) 
period 



29 
116 

58 



D 



, path = 1 • 12 • 1 • 223 • 42 • 23 • 1, configuration = ((431), (32), (2), (1)) 
Z LCM of = period 



1, 


39 


234 


234 


26 


234 


117 


117 


7 ' 


17 ' 


17 ' 


3 ' 


17 ' 


11 ' 


11 


1, 


39 


117 


117 


13 


117 


117 


117 


5 ' 


20 ' 


20 ' 


2 ' 


20 ' 


19 ' 


19 


1, 


13. 


26 


26 


26 


26 


13 


13 


7 ' 


7 ' 


5 ' 


7 ' 


3 ' 


3 


1, 


39 


468 


468 


52 


468 


117 


117 


2 ' 


71 ' 


227 ' 


11 ' 


149 ' 


31 ' 


31 



234 
117 

2G 

468 



For instance in the third example, configuration = ((43111), (4)) means that m\ 



o (1) (1) i (2) 

3, 7713 = m\ = 1, m\ 



— 1 and all the other 



,(«) 



s are 0. 



5. Size of orbit 



In this section we only consider the basic periodic A n automaton, i.e., B = Bf^. In 
addition to the period under the time evolutions, the Bethe ansatz at q = also leads to a 
formula for the size of certain orbits in the periodic automaton. Recall the quantity 



(22) 



n L (m) =detF JJ 



1 



(a) 



(a) 



» 



1 



» 



e z 



obtained in JKNj eq.(3.2) (denoted by R(v, N) therein) as the number of off-diagonal solu- 
tions to the string centre equation. Here (J) = s(s— 1) ■ • ■ (s— t+l)/t\, and L and m = (m^) 

enter the right hand side through l|ll|) and l|17|) with vy' — LS a iSji. In this special case 
the general identity known as the combinatorial completeness of the Bethe ansatz at q = 
( KN Corollary 5.6) reads 



(23) (a* 



m)x 1 x\ 



n+l' 



{la 



(a) 



(ah 



The left hand side is the character of B — Bf . The sum extending over all ; £ Z>o 
cancels out except leaving the nonzero contributions exactly when L > qi > • • • > q n . For 
example when n = 2, one has fi 6 (((3), (1))) = 6, 6 (((21), (1))) = 36, n 6 (((lll), (1))) = 18 
summing up to ( 3 j J - 60 for {q u q 2 ) = (3, 1), whereas fi 6 (((l), (3))) = 6, n 6 (((l), (21))) = 
— 18, f2 6 (((l), (111))) = 12 cancelling out for ((71,(72) = (1,3). In this sense f2z,(m) gives a 
decomposition of the multinomial coefficients according to the string pattern to. It is known 
( KN Lemma 3.7) that f2z,(m) S Z>i for any configuration, namely under the condition 
> for all (a,j) G H. Moreover it was pointed out in K OTY| that the expression ll'L'li 
for simplified by l|19() coincides exactly with eq.(2.3) in |YYT| . which is the number of 
automaton states that contain rrSp solitons with length j. Thus it is natural to ask what is 

being counted by (|22|l for the basic periodic A„^ automaton in general. 

To deal with this problem we need to consider a more general class of time evolutions. 
Let B a >i be the crystal of the Kirillov-Reshetikhin module [RMN . The crystal so far 
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written as Bi is B 1,1 in this notation. Given a pathp e5 = Bf L , seek an element v a ^ G B a ^ 
such that v aj (8) p ~ p' ® u"'-? for some p' under the isomorphism B 01 ^ ® £? ~ B ® £> aj . If 
such a t;°' J exists and p' is unique even when i;"'- 7 is not unique, we denote the p' G B by 
Tj a \p). Otherwise we say that T^ a \p) does not exist. We call p evolvable if T^ a \p) exists 

for all members of T := {T$ a) | 1 < a < n, j G Z>i}. For such a p write Tp = U a j T j a) (p)- 
We say that p is cyclic if all the paths p, Tp, T 2 p, . . . are evolvable. These paths form an 
orbit Orb(p) := Ut>o^"*(p)' which is necessarily a finite subset in B. As in the previous 
section, we let Ph(m) Q B = Bf L denote the set of highest paths whose configuration is 
m = (m^ 1 ), . . . m^™)). (Thus Ph(m) is dependent on L.) 

Conjecture 2. Given a configuration m = (m^, . . . , m^), one has two alternatives; all the 
paths in Ph(jn) are cyclic, or all the paths are not cyclic. In the former case, the following 
formula is valid: 

(24) Q L {m) = \ |J Orb(p)|. 

peP h (m) 

All the highest paths with length L < 5 are cyclic. The smallest example of non-cyclic 
Ph{m) emerges at L = 6, which is Ph(((22), (2))) only. It consists of the unique highest 

(2) 

path p = 112233, which is evolvable but not cyclic. In fact one has T{ '(p) — 213213 but in 
the next step [13] ® (213213) ~ (311223) <g> [13] whereas [23] <g> (213213) ~ (223311) ® [23]. 
Here [13] G B 2 ' 1 stands for the column tableau of depth 2, etc. Thus p' in the above 
sense is not unique, meaning that 213213 G Tp is not evolvable hence p is not cyclic. For 
L = 7, again /\(((22), (2))) is the unique case consisting of non-cyclic paths. We have 
checked the conjecture up to L = 8, where there are 5 non-cyclic ones out of 56 possible 
configurations. Some examples of Conjecture 2 are presented in the following table. (We 
write Orb(m) = U P ep h ( TO ) ( - ) rb(p), which also depends on L.) 



L 


m 


n L (m) = |Orb(m)| 


6 


((3)) 


6 


6 


((21), (1)) 


36 


6 


((1111), (11), (1)) 


12 


7 


((31), (1)) 


56 


7 


((221), (21), (1)) 


63 


7 


((21H), (21), (1)) 


133 


7 


((2111), (111), (11), (1)) 


112 


8 


((111111), (1111), (11)) 


4 


8 


((2211), (211), (11), (1)) 


192 


8 


((21111), (211), (11), (1)) 


304 



For example in the third case L — 6, m — ((1111), (11), (1)), one has 

P h (m) = {121234, 123124, 123412}, 
(25) Orb(m) = {121234, 123124, 123412, 124123, 212341, 231241, 

234121, 241231, 312412, 341212, 412123, 412312}. 



BETHE ANSATZ AT q = AND PERIODIC BOX-BALL SYSTEMS 



11 



6. Summary 

In this paper we have constructed new periodic soliton cellular automata and studied 
them by a novel application of the Bcthe ansatz. Section [21 contains the definition of the 
periodic automata associated with any non-exceptional affine Lie algebra g„. Local states 
range over the crystal B\ of Q n and the time evolution Q is a translation in the extended 
affine Weyl group. In Section we have shown that Bethe eigenvalues at q — become a 
2"P;-th root of unity, where Vi is explicitly given by the formula (|16|l . In Section 0] Vi is 
conjectured to yield the dynamical period of the An automata if F and F[b, k] in Ijltjfl are 
specified by the combinatorial Bethe ansatz jKKRl IKRj . In Section [SI the Bethe ansatz 
character formula (|23() is found to measure the size of orbits of the automata as in Conjecture 
2. 
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